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ABSTRACT 

We present VLA observations of the distribution and kinematics of the HI gas in 
the classical ring galaxy VII Zw 466 and its immediate surroundings. The H I gas corre- 
sponding to the bright optical star forming ring exhibits the typical profile of a rotating- 
expanding ring. The systemic velocity of the galaxy is found to be 14,468±25 kms -1 . 
A formal fit to the HI kinematics in the ring is consistent with both ring rotation and 
expansion (expansion velocity 32 kms -1 ). However, HI in the northeast quadrant of 
the ring is severely disturbed, showing evidence of tidal interaction. In addition, fingers 
of H I gas extend from the galaxy to the east in the general direction of the two major 
companions. We also detect a hydrogen plume from the southern edge-on companion 
galaxy (G2) pointing towards the ring galaxy. This, and other peculiarities associated 
with G2 suggest that it is the intruder galaxy which recently collided with VII Zw 466 
and formed the ring. Numerical hydrodynamic models are presented which show that 
most of the observed features can be accounted for as a result of the impact splash be- 
tween two gas disks. The resultant debris is stretched by ring wave motion in the bridge 
and later forms accretion streams onto the two galaxies. 

Finally, we detect H I emission from two previously unknown dwarf galaxies located 
northeast and southeast of VII Zw 466 respectively. This brings the total number of 
members of the VII Zw 466 group to five. Using the projected mass method, the upper 
limit of the dynamical mass of the group was estimated to be M Q = 3.5 x 10 12 Mq, which 
implies that the mass to light ratio of the group is (M/L) group w 70. This rather low 
value of M/L, as compared with other loose groups, suggests the group may be in a state 
of collapse at the present time. Plunging orbits would naturally lead to an enhanced 
probability of head-on collisions and ring galaxy formation. 



1 Visiting astronomer at NRAO. The National Radio Astronomy Observatory is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 
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1. Introduction 

The idea that galaxies might be transformed as a 
result of interactions, collisions and mergers was pre- 
sented by Toomre & Toomre (1972) to explain the 
diversity of forms seen in many of the images of pecu- 
liar galaxies (Arp 1966, Vorontsov-Velyaminov 1977). 
Observations since that time have further strength- 
ened the view that interactions play a crucial role in 
galaxy evolution. In particular, as Toomre has sug- 
gested, they may be responsible for the formation of 
elliptical galaxies through major mergers and colli- 
sions. Observations of rich clusters of galaxies at red- 
shifts of between 0.3 and 0.6 seem to show an in- 
creased fraction of blue galaxies, with a significant 
fraction undergoing an interaction (Butcher & Oem- 
ler 1984; Lavery & Henry 1988; Dressier et al. 1994). 
For those blue galaxies which are interacting, there is 
increasing evidence from HST observations that many 
have shell-like or ring-like morphologies (Lavery et 
al. 1995; Oemler 1995). Although these may not all 
be collisional systems, it is interesting to study low- 
redshift galaxies which show similar morphologies and 
ask whether these are simply analogs of processes oc- 
curring at higher redshift. In this category we would 
include the shell ellipticals (Malin & Carter 1980), the 
galaxies with ripples (Schweizer & Seitzer 1988) and 
the collisional ring galaxies (Lynds & Toomre 1976; 
see Appleton & Struck-Marcell 1996 for review ; here- 
after AS96). 

The collision of two gas-rich disk systems of un- 
equal mass ratio is probably one of the most challeng- 
ing problems in numerical modeling today (Barnes 
& Hernquist 1992, Struck 1996a). The treatment of 
shocks in the disks of two inter-penetrating galaxies 
is quite difficult to simulate and very few examples 
exist in the literature. Struck (1996b) has recently 
attempted to model both the hot and cold gaseous 
component during a head-on collision to produce the 
Cartwheel ring galaxy and many complicated physi- 
cal processes are at work during and shortly after the 
penetration of the "target" galaxy by the "intruder" . 
It is therefore of considerable interest to map the cool 
and hot components of the gas in a real collision and 
make comparisons with models. 

The ring galaxy VII Zw 466 was discussed very 
early on in the debate about the origin of these pecu- 
liar rings (Freeman & De Vaucouleurs 1974; Theys & 
Spiegel 1976; Lynds & Toomre 1976). An early pho- 
tometric study of VII Zw 466 was made by Thompson 



& Theys (1978) in which the colors and knots within 
the ring were investigated. This galaxy has also been 
the subject of a recent optical and infrared study 
(Marston & Appleton 1995; Appleton & Marston 
1995) and strong radial color gradients were found 
within the ring. 

Throughout this paper, we will adopt a Hubble 
constant of 80 kms -1 Mpc -1 . We have determined 
the heliocentric optical velocity of the galaxy to be 
14,468 kms -1 and therefore assume a distance to the 
VII Zw 466 group of 180 Mpc. 

2. Observations 

The observations were made on December 3 1994, 
using all 27 telescopes in the C configuration of the 
VLA0 . The correlator was set in a two IF (Intermedi- 
ate Frequency) mode (AC BD) with on-line hanning 
smoothing and 32 channels per IF. For each IF we 
used a bandwidth of 3.125 MHz. This provided a 
frequency separation of 97.6 kHz per channel, which 
corresponds to 22.66 kms -1 in the rest-frame of the 
galaxy using the optical definition of redshift. To 
achieve a wider velocity coverage, IF1 was centered 
at 14,131 kms" 1 and IF2 at 14,538 kms -1 . As a re- 
sult, the velocity coverage of our observations was 974 
kms -1 . A total of 4 hours 47 min was spent on source. 
Flux and phase calibration was performed using the 
sources 3C286 and 01311+678 (B1950) respectively. 

These data were first amplitude and phase cali- 
brated and bad data due to interference were flagged 
and ignored by the AIPS software. An image cube 
was created from the UV data by giving more weight 
to those baselines which sampled the UV plane more 
frequently (so-called natural weighting). This pro- 
vided a synthesized beam with a FWHM of 23.3" x 
21.9" for IF1 and 21.6" x 19.3" for IF2. 

The subtraction of the continuum emission in each 
line map was performed using a standard interpola- 
tion procedure based on four continuum maps free 
from H I at each end of the band. The resulting rms 
noise per channel was 0.27 mJybeam- 1 . The high- 
est dynamic range in any channel map was 5:1 (at 
v=14,402.8 km/s) with a peak flux of 1.35 mJy/beam. 

In order to determine the total H I distribution we 
used the following technique. Initially, we smoothed 

2 The Very Large Array is operated by Associated Universities 
Inc. under cooperative agreement with the National Science 
Foundation. 
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all channel maps to a resolution 44.0" x 44.0", twice 
that of the synthesized beam. New maps were formed 
comparing, pixel for pixel, the original full resolu- 
tion maps to the smoothed ones. The pixel values 
of original maps were copied to the new ones only 
if the signal-to-noise ratio of the smoothed map at 
that point exceeded 2. The total HI surface density 
map was produced by adding the new maps together. 
The same technique was used to create the first and 
second moment maps of the distribution. Using this 
procedure we effectively give more weight to points 
associated with low surface brightness emission. 

Optical observations of a 16x16 arcmin field were 
obtained during photometric conditions on April 2 
1995, with the Fick Observatory 0.6 m telescope and 
CCD system (Applcton, Kawaler & Eitter 1993). A 
500 sec exposure using an R-band filter was taken. 
Calibration of the photometry was performed using 
data taken from the star G10 50 (Landolt 1992) 

3. The VII Zw 466 Group and its H I Distri- 
bution 

In Figures 1 and 2 we show the VII Zw 466 group 
and its environment. The ring has two nearby com- 
panions, labeled Gl and G2 as well as a background 
galaxy Bl (Appleton & Marston 1995). Galaxy Bl 
(first noted as background by R. Lynds; see Theys & 
Spiegel 1976) has a velocity of V Q = 25,042 kms" 1 , 
determined from our own unpublished optical spectra. 
Also marked in Figure 1 are two new dwarf galaxies, 
G3 and G4, discovered during the H I mapping of the 
group, presented in this paper. 

In Figure 2 we show the integrated neutral hydro- 
gen emission from the VII Zw 466 field superimposed 
on a B-band image of the group (from Appleton & 
Marston 1995). The emission from the ring galaxy is 
highly disturbed. The brightest emission comes from 
the northern quadrant of the ring and sweeps around 
in a clockwise direction along the western quadrant of 
the ring. A filament of H I extends from the ring in a 
south-easterly direction, pointing to a region between 
Gl and G2. In addition, there is a marked depression 
in the H I emission from the south-eastern quadrant of 
the ring. It is in this quadrant that the optical emis- 
sion in the ring is rather faint. The H I distribution 
in VII Zw466 is similar to that seen by Higdon (1993; 
1996) in the Cartwheel ring galaxy, with scattered HI 
extending in the direction of the three possible com- 
panions. As with VII Zw 466, the Cartwheel plumes 



are highly asymmetric, with all the gas outside the 
galaxies being found on one side of the ring, highly 
suggestive of a hydrodynamic splash resulting from 
the gas on gas collision between the galaxies. 

H I emission is also detected from the edge-on com- 
panion G2. The HI has a large peak centered on 
the south-eastern end of the major axis with fainter 
emission extending along the major axis. A very faint 
filament extends from G2 back towards the ring. 

No emission down to a level of 0.68 mJybeam -1 
(2.5c) was detected from the peculiar elliptical galaxy 
Gl. If we adopt a typical velocity width of 300 km s -1 
for the assumed velocity of any gas in the elliptical 
galaxy, this limit translates into a 2.5a upper limit of 
3.1 x 10 9 Mq for the hydrogen mass of the companion. 

We also detected gas from two new group members 
G3 and G4 located 3.2' and 4.2' to the southwest and 
northwest of VII Zw466, respectively. At the distance 
of the group these correspond to projected linear sep- 
arations of 166 kpc and 220 kpc from the ring. 

4. Integrated Properties 

In Table 1 we present the HI and optical proper- 
ties of VII Zw 466, its two major companions (Gl and 
G2), and the two newly discovered dwarf galaxies, G3 
and G4. The R-band luminosites, Lr, of Gl and G2 
are from Appleton & Marston (1995) while the op- 
tical properties of G3 and G4 were determined from 
the Fick CCD images. There is excellent agreement 
between the photometry of the galaxies in common 
between the previously published KPNO magnitudes 
and the Fick magnitudes. 

The velocity profiles of the neutral hydrogen dis- 
tribution of VII Zw 466, G2 and two dwarfs are pre- 
sented in the Figures 3-6. The H I profile of VII Zw 466 
presents the typical two-horned profile of a rotating 
disk. We estimate the systemic velocity of the ring 
to be 14,468 ± 25 km sec -1 , in close agreement with 
the optically determined velocity of 14,490 km sec" 1 
(Theys & Spiegel 1976). However, the high velocity 
horn is not completely symmetrical. Emission seen at 
velocities higher than 14,520 km sec -1 is associated 
with the gas plumes outside of the visible galaxy. 

The HI profile of G2 is single-peaked with a AVx/ 2 
= 84 km sec -1 . The galaxy appears to be edge-on, 
so for all calculations we used an inclination «=90°. 
The two dwarfs, though, have more peculiar HI pro- 
files. The velocity profile of G3 exhibits a low-level 
asymmetry around 14100 km sec -1 which may not be 
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real. For this reason, we did not attempt to calculate 
AV1/5 f° r this galaxy. Finally, G4, the least massive 
galaxy in the group, has a very narrow H I profile with 
AV 1/2 = 25.3 kmsec^ 1 . For both G3 and G4 we as- 
sumed an inclination of i=45°, since the disks of the 
galaxies were only marginally resolved in the optical 
images of the group. 

We also include in Table 1 a rough estimate the 
dynamical mass of each group member. The val- 
ues of Md corresponding to the galaxies G2, G3 and 
G4 were calculated using the formula: 

(hAVi cosec(i)) 2 R HI 
M d = ^ g — 

which assumes that the gas is in bound circular orbits. 
Values of the H I radius were obtained by inspection 
of the run of emission centroid with radius for each 
galaxy and in all cases the H I disks were resolved and 
a radius determined. 

The estimate of the dynamical mass of VII Zw 466 
was produced by the fit to the rotation of the HI 
ring discussed in Section 6. For completeness, the 
dynamical mass of Gl was simply estimated from its 
observed luminosity, using an R-band mass to light 
ratio of M/L=15, typical of elliptical galaxies (Faber 
& Jackson 1976). Finally, we present the mass to 
light ratio M d /L R of G2, G3, G4, and VII Zw 466. 
We note the surprising small values of M d /L R as well 
as Mhi/Lr for the spiral galaxy G2. We believe that 
this may be related to the tidal stripping of its mass. 
We will elaborate more on that in Sections 6 and 8. 

5. The Dynamical Mass of the VII Zw 466 
Group 

Our observations indicate that there are now four 
galaxies G1-G4 which form a loose group around 
VII Zw 466. Since we know the recession velocities of 
those galaxies and their angular separations, we can 
use a dynamical method to get an independent esti- 
mate on the total mass of the group and subsequently, 
its dark matter content. This is the first time that the 
large-scale matter content of a group containing a ring 
galaxy has been estimated. 

The primary method used is the projected mass 
estimator proposed by Bahcall & Tremaine (1981). 
However, for comparison we also estimated the mass 
of the group using a mass weighted virial theorem 
approach. Both give similar results. 



The projected mass method is based on the idea 
that, for a dynamically bound system with one mas- 
sive object and several smaller companions, one can 
examine two extreme cases for the possible orbital 
types of the group members. The first case assumes 
an isotropic distribution in the velocities of the mem- 
bers which implies that the average value for the ec- 
centricity of the orbits is < e 2 >= 1/2. The other 
extreme is that all orbits of the group members are 
radial implying that < e 2 >— 1. These lead us to de- 
fine two specialized estimators of the dynamical mass: 

16 N 

i=l 

and 

i—l 

where v Z i is the difference in the recessional ve- 
locity between the group member i and the central 
massive object and Ri is its projected distance from 
the central object (see Bahcall & Tremaine 1981, for 
details). 

When one has no specific information on the dis- 
tribution of the eccentricities of the group one can use 
a dynamical mass estimator which is the average of 
the two extremes, that is: 

1 24 A 

Mo = ^ (Mi + M R ) = — v 2 zi Ri 

i=l 

In the VII Zw466 group the most massive object is 
the elliptical galaxy Gl. However, since its velocity 
is poorly determined (no HI detection) and its mass 
is inferred from its luminosity and an assumed mass 
to light ratio, we decided to to use the barycentric 
velocity and position of the group as the point of ref- 
erence, rather than Gl. Using the figures presented 
in Table 1, we find that the barycenter is located at 
a cm =12 h 32 m 11.2! and <5 cm =+66 23'58.51" (J2000) 
and its systemic velocity is w cm =14,143 kms -1 . The 
new position for the center lies only 1 galactic radius 
away from the elliptical (the elliptical still dominates). 
Our results for the mass of the group are presented in 
Table 2. 

Although the use of the barycenter rather than a 
single galaxy as the center may seem to contradict the 
assumption that the group is dominated by a single 
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central massive object, we feel that this is a com- 
promise position, given the poorly determined opti- 
cal velocity of Gl. If the group is dominated by dark 
matter, then the barycenter of the group might be 
more appropriate anyway. We note for comparison 
that had we used Gl as the center, our estimate of 
the dynamical mass would only be greater by a factor 
of 1.83 than the result presented in Table 2. 

The estimate of the total dynamical mass of the 
group is M = 3.5 x 10 12 M . The sum of the dy- 
namical masses of each individual (ith) galaxy from 
Table 1 is: 



M T = V M dl = 3.55 x 10 n M Q 

which implies that the amount of dark matter dis- 
tributed outside the members of the group is: 



M 
M T 



9.9 



Table 1 also shows that the total luminosity of the 
group is L^t=5.02 x 10 10 Lq, implying that the total 



mass to light ratio for the group is (M/L) s 



70. 



For completeness we also calculated the mass of 
the group using a method based on the mass- weighted 
virial theorem. In this approach the mass of the group 
is given by the formula: 



Mvt = 



3- 



E M < 



diV z 



2G 5 



E E 

i=l j=l,j>i 



M di M dl 



E M di 



where M d i is again the dynamical mass of a group 
member i and is the the projected distance be- 
tween the members i and j. This results to a group 
mass M VT = 6.2 10 12 M Q , which is only 1.78 times 
larger than the projected mass method estimate. 

We believe that the projected mass method esti- 
mate is more accurate than the one resulting from the 
virial theorem technique. Bahcall & Tremaine (1981) 
have shown through numerical simulations, that the 
virial theorem technique may be biased and inefficient 
in particular for systems with small number of mem- 
bers. In addition, the projected mass method has the 
added advantage that it explicitly takes into account 



the possibility that some or all of the orbits in the 
group are radial. The latter reason is very relevant in 
this study. Unlike most group studies, we know that 
in the VII Zw 466 group at least two of the galaxies 
(perhaps G2 and VII Zw 466) must have collided ra- 
dially in order to produce the ring. 

The projected mass method and the virial mass 
method are valid only if the group is gravitationally 
bound. If some of the galaxies are leaving the group, 
or are chance projections of foreground or background 
galaxies, or the group is not stable, then the resulting 
dynamical mass and group crossing times would be 
incorrect. 

Other studies of loose groups of galaxies (e.g. 
Ramclla, Geller & Huchra 1989) show that they have 
a typical (M/L) ~ 150-300, while in compact groups 
(M/L)« 50 (Hickson et al. 1992). In some scenarios, 
it is believed that compact groups are believed to be 
continously forming at the center of rich loose groups 
and are in a state of collapse. The fact that the mass 
to light ratio of the VII Zw 466 group is similar to that 
of compact groups and also the presence of a least one 
observed collisional ring galaxy pair in the group, sug- 
gests that, like the compact groups, VII Zw466 group 
may also be collapsing. 

6. The H I Kinematics of the Ring 

In order to understand the kinematical behavior of 
the H I gas in the ring of VII Zw466, we used a model 
which assumes that the ring may be both rotating and 
expanding. This model has been used successfully in 
other ring galaxies, such as the Cartwheel (Fosbury 
and Hawarden 1977; Higdon 1993). Using the optical 
B-band image of the group and assuming an intrin- 
sically circular ringf], we estimated the inclination of 
the ring to be i = 37° and the position angle of the 
major axis to be -5° (N through E). Then, we mea- 
sured the velocities of the HI gas along azimuth of 
the ring and performed a three parameter fit using 
the following function: 



v = <2q + aisin((p + <po) 



(1) 



3 We caution that we have assumed that the ring is an intrinsi- 
cally circular structure and that its elliptical shape is entirely 
a result of an inclined viewing angle. We acknowledge that 
in models of off-center ring formation, in the early expansion 
stages, the ring is somewhat non-circular. This may lead to an 
overestimation of the inclination of the ring and would require 
a more complicated treatement of the kinematic data. 
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where ao,ai and <j)Q were the free parameters. A 
rotating and expanding ring would exhibit a simple 
sinusoidal shape, with a phase offset from the ma- 
jor axis that is related to the amplitude of the ex- 
pansion. The formal fit for the ring resulted in a 
X 2 =5.61. However, the points of the northeast quad- 
rant of the ring exhibit a peculiarity related to one of 
the HI plumes, and the errors associated with them 
are larger. If we exclude those points, (which is rea- 
sonable since we are only trying to model the ring) 
the fit improves considerably, giving a x 2 =2.77 for 
a =14,500kms- 1 , ai=83kms -1 , and0 o =274°. This 
fit is presented in Figure 7. The azimuthal component 
of the velocity at the ring radius and the radial expan- 
sion velocity, corrected for inclination are found to be 
v ma a; = 137 kms -1 and w exp =32 kms -1 . We note that 
this expansion velocity is significantly larger than the 
8 kms -1 obtained optically by Jeske (1986). We be- 
lieve that the combination of the two-dimensional na- 
ture of our fit and the high velocity resolution of our 
observations makes our measurement more reliable. 

In Figure 8 we also present the channel maps of the 
HI observations. In the rest frame of the galaxy, the 
channel separation is 22.66 kms -1 . HI emission from 
the VII Zw 466 is detected for the first time around 
14,561 kms -1 and originates from the southern side 
of the galaxy. As we move to lower velocities, the 
emission is typical of a rotating ring. We observe 
that the emission splits in two centroids located to 
the east and west of the optical ring. However, the 
eastern centroid does not overlap spatially with the 
optical picture of the ring. It resides at the outside 
of the optical ring structure and at 14,470 kms -1 is 
deformed into a small plume (hereafter "Plume A") 
pointing to the southeast. The amount of HI gas 
associated with this plume is M^=1.6 x 10 8 Mq. 
The emission from the western centroid does coincide 
with the optical ring. It is also considerably stronger, 
indicating large quantities of gas. Strong emission 
from the west side of the ring was also detected in the 
radio continuum observations of the galaxy (Ghigo & 
Applcton, in preparation). 

As we proceed to channel maps with velocities be- 
low 14,470 kms -1 , the behavior of the gas is very dis- 
turbed. Firstly, we notice that at 14,448 kms -1 there 
is practically no H I gas associated with the ring. This 
is very peculiar since there is a strong emission in the 
preceding and following channel maps. Also, there is 
a filament (hereafter "Plume B" ) of H I gas present 
in three consecutive channel maps (14,402 kms -1 - 



14,357 kms -1 ). This is the strongest emission that we 
observe. The filament initially points to the northeast 
of the ring and as we progress to lower velocities, it 
points to the east-southeast of the ring in the gen- 
eral direction of the companions Gl and G2. We cal- 
culated that the HI gas associated with plume B is 
M B =4.2 x 10 8 M . 

At this point, it also useful to discuss the HI emis- 
sion from the companion G2. The emission from 
G2 is located in the channel maps between 14,516- 
14,334 kms -1 . However, we notice HI gas at a veloc- 
ity of 14,584 kms -1 and a lack of HI gas at 14,425 
kms -1 . Moreover, the emission centroid from G2 in 
the channel maps after 14,425 kms -1 is rather elon- 
gated pointing clearly towards the ring galaxy. Ta- 
ble 1 shows that the value of Mhi /Ljj for G2 is also 
unusually low for a late-type spiral, suggesting that 
it has been tidally stripped. (We also note that the 
total mass to light ratio of G2 is extremely low, im- 
plying that the dark matter halo of G2 has also been 
stripped). 

The dynamical behavior of the HI gas is consis- 
tent with G2 being the galaxy that disturbed the 
progenitor H I disk of VII Zw 466 and caused the for- 
mation of an expanding star forming ring. Numerical 
simulations of the gas dynamics of a nearly head-on 
collision between a large gaseous disk galaxy and a 
smaller companion disk galaxy (Struck 1996a;1996b; 
see also Section 8) show that large quantities of gas 
are stripped away from the disk of the target galaxy 
by the companion during the collision. Most of the 
ejected material either falls back to the large dis- 
turbed disk or it spirals towards the smaller com- 
panion which eventually accretes it. We will argue 
that the asymmetric nature of the HI distribution 
around G2 is strongly suggestive of the onset of ac- 
cretion from one such tidal stream. 

Plumes A and B and emission which extends from 
the southern end of VIIZw466 (see Figure 2) prob- 
ably form part of a two-horned HI bridge between 
VII Zw 466 to G2. This rather special morphology is 
predicted by the models and strongly suggests that 
the plumes are splattered debris from the collision. 
The total amount of HI gas in plumes A and B is 
Ma+M b =5.8 x 10 8 Mq, which is approximately 
16% of the mass of the H I disk of VII Zw 466. A more 
detailed discussion of the formation of these struc- 
tures is given in the next section. 

To further estimate the energy needed to gener- 
ate the two HI plumes, we performed an order of 



6 



magnitude calculation of the kinetic energy of the gas 
in the plumes. We used the HI masses given above 
and the relative velocities of plumes A and B Ava 
= 14 kins -1 and Avb = 88 kms -1 with respect to 
VII Zw 466. The kinetic energy associated with their 
gas was found to be KE^= 3.1 xlO 53 ergs and KE^= 
3.2 xlO 55 ergs respectively. Although projection ef- 
fects are not included, this implies the energy required 
to displace the plumes from VII Zw 466 is of the order 
of 10 55 ergs. This is approximately 10% the gravita- 
tional potential energy of the VII Zw 466 - G2 pair, 
indicating that the collisional scenario is plausible. 

7. The Possible Role of the Elliptical Com- 
panion Gl 

In our discussion so far we presented evidence in fa- 
vor of the edge-on disk- like companion G2 as the most 
likely intruder galaxy. Can we completely rule out 
the possibility that Gl, the brighter elliptical galaxy, 
is the galaxy that plunged through the center of the 
progenitor disk of VII Zw466? After all, an inspection 
of Figure 1 and 2 shows that the elliptical system lies 
closer to the minor axis of the ring, is slightly closer 
to the ring (2.5 rather than 3 ring diameters away) 
and has been shown to have peculiar optical isophotes 
(Appleton and Marston 1995). In addition, its optical 
velocity (though poorly known) differs from the ring 
velocity by over 350 kms -1 , suggesting a high impact 
velocity with respect to the ring. 

Although we cannot completely rule out the possi- 
bility that the elliptical galaxy is the intruder galaxy, 
we believe the evidence at present is against it. Firstly, 
all our simulations to date lead to the formation of an 
accretion stream, which is a bridge between the two 
galaxies involved in the collision. In the only similar 
21cm HI study so far published, the same behavior 
is found for the Cartwheel galaxy, where a long HI 
streamer connects the Cartwheel to its most distant 
disk companion (Higdon 1996). The fact that the el- 
liptical system is not connected to the ring, nor shows 
any evidence for HI down to a low column density (see 
Table 1) seems to rule out accretion of cool gas onto 
galaxy Gl. A careful analysis of the color map of the 
VIIZw466 system (see Appleton and Struck-Marcell 
1996) seems to rule out any evidence for a dusty disk 
inside Gl, similar to those seen in other ellipticals in 
which accretion is believed to be taking place (e.g. 
Sparks et al. 1985). Is it possible that any accreting 
gas might have become ionized? If so, such emission 



must be weak, since imaging in Ha by Marston and 
Appleton (1995) failed to show any evidence for ion- 
ized hydrogen. In addition, unpublished radio contin- 
uum maps of the system by Ghigo and Appleton (in 
preparation) show that both VII ZW 466 and G2 are 
strong radio emitting objects at a wavelength of 6cm, 
but that no emission is observed from the elliptical 
galaxy Gl. Similarly, Thompson and Theys (1978) 
found comparably blue colors in VII Zw 466 and G2, 
but not in Gl. These points again suggest that Gl is 
a rather inactive object and probably has not suffered 
accretion recently. 

Perhaps the strongest arguments in favor of Gl be- 
ing the intruder relates to the position of the galaxy 
relative to the minor axis of the ring and its high rela- 
tive velocity compared with the ring^. The galaxy Gl 
is only a little closer to the ring than G2, and so from 
this point of view neither is favored. N-body simula- 
tions by Huang and Stewart (1988) and Appleton and 
James (1990) show that very respectable rings can be 
generated in collisions of quite low inclinations be- 
tween the intruder and the target disk, reducing the 
need to always find the intruder close to the minor 
axis of the ring in all cases. Even the high relative 
velocity is not a strong argument for identifying Gl 
as the intruder. Again N-body simulations of ring 
galaxies (e. g. R. A. James, private communication) 
show that moderately unbound collisions are turned 
into bound systems via dynamical friction and that 
it is likely that the relative velocity of the intruder 
would be significantly reduced after it has travelled 
a few ring diameters from the target (as observed). 
In conclusion, we believe that the HI data provide 
strong evidence that the disk galaxy G2, rather than 
the elliptical Gl, is the likely intruder in this system. 

8. Models 

The HI distribution of Figure 2 suggests that we 
are witnessing the messy aftermath of a collisional 
splash between two galaxies. We present here a nu- 
merical simulation which supports this picture, and 
provides further details about the dynamics. The 
model presented here is not a detailed model of this 
particular system, but should be viewed as a working 
example of how the collisional scenario (disk plus disk 

4 Indeed, although Theys and Thompson (1978) seemed to favor 
G2 as the intruder based on color evidence, they could not 
rebut this argument in a convincing way (They appealed to 
three-body effects). 
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collision) can produce much of the observed struc- 
tures. Until higher resolution HI observations are 
available, a more detailed model which tries to fit 
all aspects of the system perfectly is not yet justified. 
The model presented here is part of larger program of 
simulations between two gaseous galaxy disks. Some 
preliminary results have been presented in AS96 and 
Struck (1996a), and a journal paper on these simula- 
tions is in preparation (Struck 1996b). 

The simulation code is described in the above ref- 
erences, and an earlier version is discussed in Struck- 
Marcell & Higdon (1993). In this code, the large-scale 
collision dynamics are modeled within the restricted 
three-body approximation, which is adequate for the 
early transient stages of galaxy collisions (e.g. Gerber 
& Lamb 1994). In these models each galaxy consists 
of a rigid halo and a gas disk. The potential of the 
primary galaxy is such as to give a rising rotation 
curve, which becomes nearly flat at large radii, and 
extends to several disk diameters. A simple softened, 
point-mass form is used for the companion potential. 
The gas dynamics is computed with a smooth particle 
hydrodynamics (SPH) algorithm with a spline kernel. 
In addition to the force from the large-scale potential, 
a local self-gravity is computed between neighboring 
gas particles within the disks, which allows the forma- 
tion of clumps or clouds. Shear dominates over local 
gravity on scales larger than that on which the latter 
is calculated. 

Simple models for heating and cooling are also in- 
cluded in this simulation. However, isothermal com- 
parison runs show that the large-scale morphology is 
little effected by these terms (see Struck 1996b for de- 
tails). The most important result of including these 
terms, for present purposes, is that the simulations 
show a prompt, adiabatic cooling of bridge material. 

The two disks were initialized in centrifugal bal- 
ance, with the companion set several primary-disk di- 
ameters from the target, and with a slight sideways 
velocity to allow a slightly off-center impact. The 
x-y plane is the initial midplane of the primary disk, 
while the companion midplane lies on a translated y-z 
plane. The companion orbits in the x-z plane. There 
is sufficient time before impact for the heating/cooling 
terms of both disks to reach a quasi-equilibrium. 

At the time of impact, extremely strong shocks 
form. A good deal of primary material is splashed 
out, so that the bridge ultimately contains a fairly 
even mix from both progenitors. Similarly, much gas 
from the companion is left behind. Although this 



simulation was not specifically produced to model 
VII Zw 466, the model and timestep shown in Fig- 
ure 9 were chosen as a reasonably good match. In 
a comparison model using a gasless companion but 
otherwise identical parameters to the current model 
show that no gas bridge was formed and little gas 
debris was ejected from the gas disk. Thus, we see 
that collisional hydrodynamics (especially a gas on 
gas collision) seems to plays a key role in the forma- 
tion of the double-horned HI bridge, in contrast to 
bridges formed in prograde, planar interactions (e.g. 
Toomre & Toomre 1972). VII Zw 466 is apparently 
an interesting example of such a highly dissipative 
gas-dynamical collision. 

Figure 9a,b gives two orthogonal views of the sys- 
tem some time after impact, when the companion has 
moved a good distance away, with a prominent bridge 
stretching between it and the primary. We have not 
attempted to find an optimal viewing angle for the 
model, though it appears that one about half-way be- 
tween the two orthogonal views would be quite good. 
In any case the morphological similarities between 
Figures 2 and 9 are quite strong. First of all, both 
have a similar double-horned bimodal appearance, or 
alternately, a central hole in the bridge material. An- 
imated views of a sequence of model timesteps pro- 
vide a simple explanation. The ring galaxy itself is 
produced by first an inward gravitational impulse fol- 
lowed by a centrifugal rebound. The same kinematics 
occurs in the bridge material, with the addition of 
stretching between the galaxies. Thus, in this picture 
the bridge is essentially a ring stretched the vertical 
direction, or a cylindrical annulus. 

Figure 9a,b also shows the pattern of infall onto the 
two galaxies. In the case of the primary disk there is 
significant infall out of both streams, and the accre- 
tion tends to be concentrated in the central regions. 
In the case of the companion the (highly supersonic) 
accretion out of one stream is much stronger than the 
other, though the details are quite time-dependent. In 
both cases the accretion streams swirl around the nu- 
cleus and then merge into the ambient rotational flow. 
Figure 9a shows that gas swirls out a relatively large 
distance on the side opposite the main stream, which 
may account for the HI asymmetry of G2 in Figure 2. 
The models indicate that the accreting material is 
shocked after circling nearly all the way around the 
companion, and there is much hot and warm gas con- 
centrated on the main stream side. The companion is 
a radio continuum source, and the models predict that 
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this emission would be offset from the galaxy center 
if arises from such shocks. This is in fact observed 
(Ghigo and Appleton, in preparation). In the model, 
the companion disk is highly disrupted in the colli- 
sion, and largely reforms in the accretion process. Al- 
though we have no direct evidence for accretion onto 
VII Zw 466 itself, we note that the large circular loop 
seen at optical wavelengths in the north-west quad- 
rant of the ring might be the result of bridge material 
falling back onto the disk. Such material might induce 
star formation in a ring, which is consistent with Ha 
observations of the loop (Marston & Appleton 1995). 

In an offset collision, the primary disk receives an 
azimuthal swing impulse as well as the radial one. 
This leads to asymmetric and offset spatial distribu- 
tions as seen in both Figure 2 and 9b. Figure 9b 
also shows a developing oval ring rather like that in 
VII Zw 466. This is deceptive, however, since it is the 
second ring in the model, with the first ring moving 
out of the disk. At present there is no independent 
evidence to suggest that the ring in VII Zw 466 is not 
the first. Thus, it appears that the impact velocity 
in VII Zw 466 was high relative to orbital or epicyclic 
speeds in the primary disk, because the companion 
has moved several ring diameters. The model shows a 
similar morphology with the first ring, but the bridge 
has hardly developed at that time. With the aid of 
kinematical constraints from higher-resolution obser- 
vations at 21cm, future models for this system should 
be able to resolve this timing problem, and provide 
more specific predictions. 

9. Conclusions 

In the previous sections we presented our H I obser- 
vations of the VII Zw 466 group. Though only at mod- 
erate resolution due to the large redshift (z=0.048) of 
the group, we were able to draw the following conclu- 
sions : 

• The H I emission from VII Zw 466 is consistent 
with that of a rotating and expanding ring. 
A simple kinematical model fits the ring suffi- 
ciently well and indicates a rotation curve with 
v max — 137 kms" 1 at the radius of the ring, and 
an expansion velocity v exp =32 kms -1 . How- 
ever, one quadrant of the ring has anomalous 
velocities and is associated with disturbed HI 
outside the ring. 

• The northern part of the ring shows evidence of 



the tidal interaction. Plumes of HI gas (5.8 x 
10 8 M , or 16% of the total H I disk) lie outside 
the ring and point in the general direction of 
the two nearby group members. The dynamics 
of this gas is consistent with a collisional origin 
for the ring. 

• Of the two nearby companions, only the edge- 
on spiral, G2, contains neutral hydrogen above 
our detection limits. Its HI properties are ab- 
normal. Not only does G2 exhibit a plume 
which points towards the VII Zw 466, but it also 
has very low values of M#//L# and M^/L^. 
This strongly suggests that G2 is the intruder 
galaxy that collided head-on with VII Zw 466 
and created the ring (Lynds & Toomre 1976). 
We conclude that the elliptical galaxy Gl has 
played little role in the ring formation, despite 
its unusual "boxy" optical isophotes (Appleton 
& Marston 1995). 

• Models suggest that the plumes represent a mix- 
ture of gas ejected from the progenitor of the 
ring galaxy and material stripped from the sus- 
pected intruder galaxy. The hydrodynamical 
models provide a good match to the plume mor- 
phology and suggest that accretion onto both 
galaxies is already underway. The models also 
suggest the possibility that the companion disk 
was strongly disturbed in the collision, and that 
the observed asymmetries are the result of par- 
tial reformation of the (G2) disk in the accretion 
process. 

• We identified two new members of the extended 
VII Zw 466 group. These galaxies are dwarfs 
and are detected for the first time in H I. 

• An upper limit to the dynamical mass of the 
group was estimated using two different tech- 
niques and was found to be in the range M a 
= 3.5 to 6.2 x 10 12 M . Adopting the lower 
end of the range as being more reliable, this im- 
plies that the dark matter not associated with 
the group members is 10 times greater than the 
sum of the dynamical masses of all group mem- 
bers. The resulting mass to light ratio for the 
group is (M/L) group « 70. 

• The estimated M/L for the group is roughly a 
factor of three times lower than the median M/L 
for loose groups, but is similar to that estimated 
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for compact groups of galaxies. As with com- 
pact (Hickson) groups, the low M/L ratio may 
imply that the group is not in equilibrium, but 
may be collapsing. Such a scenario would natu- 
rally give rise to an increased incidence of deeply 
penetrating orbits and might increase the prob- 
ability of a head-on collision leading to the for- 
mation of a ring galaxy. 

We note that, if the numerical models are correct, 
the gas plumes seen in these H I observations may rep- 
resent only a part of the debris formed as the intruder 
collided with the disk of VII Zw 466. If the collision 
was quite recent (as implied by the fact that only one 
ring is seen in VII Zw 466 - see AS96 for review) , then 
it is possible that much of the shocked gas accreting 
onto the two galaxies is still strongly ionized. Deep 
Ha observations or high sensitivity UV or X-ray ob- 
servations may be fruitful. 
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Table 1 

Properties of the VII Zw466 Group 





VIIZw466 


Elliptical (Gl) 


Spiral (G2) 


Dwarf- 1 (G3) 


Dwarf-2 (G4) 


a(J2000) 


12 h 32 m 05.27: 


1 2^32^13.30! 


12h32 m 11.95? 


12 h 31 m 37.99! 


12^31^35.11! 


c5(J2000) 


+66°24'13.65" 


+66°23'59.80" 


+66°23'19.92" 


+66°22'35.56" 


+66°26'55.26" 


Vhi (kms- 1 ) 


14,468 ± 25 


[14,100] a 


14,457 ± 25 


14,046 ± 25 


14,223 ± 25 


AV 1/2 (kms- 1 ) 


202 




84 


113 


25 


AV 1/5 (kms- 1 ) 


256 




119 




40 


Rh i (arcsec) 


11 




6.7 


9.3 


4.6 


/ S(v)dv b 


0.541 


< 0.408 c 


0.155 


0.383 


0.092 


M m /M Q d 


4.1 x 10 9 


< 3.1 x 10 9c 


1.2 x 10 9 


2.9 x 10 9 


7.0 x 10 8 


R-band mag. e 


14.94 


14.64 


15.62 


16.18 ± 0.05 


17.67 ± 0.10 


M/ (M ) 


4.10 xlO 10 


3.00 xlO 11 


2.31 xlO 9 


1.16 xl0 1Og 


1.80 xl0 8g 


L* (L©) 


1.6xl0 10 


2.0 xlO 10 


8.2 xlO 9 


4.9 xlO 9 


1.3 xlO 9 


M HI /L R (M /L Q ) 


0.26 


< 0.16 


0.14 


0.59 


0.56 


M d /Lfl (M /L ) 


2.59 


15 


0.28 


2.36 


0.21 



a optical velocity from Thcys & Spiegel (1976) 

b where J S(v)dv is the integrated flux density in Jy kms -1 

c assumes a hypothetical velocity width of 300 kms -1 and a solid angle equal to two VLA beams 
d where M HI /M Q = D 2 x 2.356x 10 5 / S(v)dv and D = 180 Mpc 
e from Appleton & Marston (1995) 

(0. 5AVi cosec(i)) 2 R H T 

t Md = — q is an estimate of the total dynamical mass of the galaxy. 

s since the dwarf was unresolved from the B-band image we assumed an inclination i = 45° 
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Table 2 

Parameters for the Dynamical Mass for the VII Zw 466 Group 



Galaxy R t a (kpc) t^kmsj g, = (M Q ) 



VII Zw466 


33.7 


325 


8.00 x 10 11 


Gl 


11.0 


-43 


4.57 x 10 9 


G2 


33.9 


314 


7.51 x 10 11 


G3 


188.4 


-97 


3.98 x 10 11 


G4 


244.1 


80 


3.51 x 10 11 



a Projected distance from the barycenter 

b Velocity offset from the systemic velocity of the barycen- 
ter (w cm =14,143 kms -1 ). 
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Fig. 1. — Grey-scale image of VII Zw 466 through a V-band filter taken at Fick Observatory. The positions of the 
two newly discovered dwarfs are marked as G3 and G4. 



Fig. 2. — Grey-scale image of VII Zw 466 through a B-band filter overlayed with a contour map of the integrated 
HI distribution. The contour increment is 57.15 Jybeam -1 ms" 1 and the level of the first contour is also 57.15 
Jybeam -1 ms -1 . 



Fig. 3.— The global HI profile of VII Zw 466. 



Fig. 4. — The global HI profile of the Companion Galaxy G2. 



Fig. 5.— The global HI profile of the Dwarf Galaxy G3. 



Fig. 6.— The global H I profile of the Dwarf Galaxy G4. 



Fig. 7. — The HI velocity as a function of azimuth around the ring. The solid circles were excluded from the fit. 



Fig. 8. — Contour plots of the 27 channel maps of VII Zw 466. The velocity of each channel is displayed in the upper 
right corner. The contour increment is 1.35 xlO -4 Jybeam -1 (0.5 a level) and the lowest contour displayed is at 
2.5cr. The solid crosses indicate the position of the ring, Gl, and G2, while the open cross indicates the position of 
the background galaxy Bl. 



Fig. 9. — Two orthogonal views of a single timestep in a numerical hydrodynamical model of the collision between 
two galaxies with gas disks, a) x-y plane, the initial midplane of the primary disk; companion is seen edge-on at 
right, b) y-z plane, the initial plane of the companion, now seen at the bottom. Both frames optimized to illustrate 
the accretion streams. Only 20% of the 19,640 particles in the primary disk, and 50% of the 3240 particles in the 
companion are plotted. Line segments proportional to the projected direction and magnitude of the velocity of 
selected particles, with the particle plotted at the base, reveal the flow. The selected particles include all particles 
originating in the companion galaxy, with temperatures of less than 10 times the initial temperature (i.e., < 30,000- 
50,000K), and located more than 0.5 (dimcnsionless) units from the z=0 plane. See references in text for further 
details. 
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